Polarized networks, diameter, and synchronizability of networks 
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Previous research claimed or disclaimed the role of a small diameter in the synchronization of a 
network of coupled dynamical systems. We investigate this connection and show that it is two folds. 
We first construct two classes of networks, the polarized networks and the random networks with a 
fixed diameter, which exhibit very different synchronizability. This shows that the diameter itself is 
insufficient to determine the synchronizability of networks. Secondly, we derive analytic estimates 
on the synchronizability of networks in terms of the diameter, and find that a larger size of network 
admits of a more ffexible synchronizability. The analysis is confirmed by numerical results. 
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I. INTRODUCTION 



Since some critical properties of complex networks 
were revealed [l|, both theoretical and experimental in- 
vestigations of complex networks have received rapidly 
increasing interests |2|- After earlier verification of cer- 
tain characteristics in various computer-generated and 
real- world networks Q], the focus of research shifted to 
the dynamics on complex networks [Jl- The synchro- 
nization of a network of dynamical systems has been ex- 
tensively analyzed and has found many applications in 
secure communication, distance detection, and other en- 
gineering problems |^ . 

It is now well known that the synchronizability of a 
network can be characterized by the Laplacian spectrum 
of its associated graph. Consider a network of coupled 
dynamical systems x* = F(x*)-|-(T^"^j LijH(x-'), where 
X* are the variables associated with node z; F and H are 
evolution and output functions, respectively; a is control 
strength; and Ly are entries of the Laplacian matrix, 
defined by La = ki, the degree of node i, Ly = — 1 
if nodes i and j are connected, and Lij = other- 
wise. Order the eigenvalues of the Laplacian matrix as 



= Ai < Aa < 



< A„. Then, the larger the two 



quantities A2 and Aa/A^ are, the better synchronizabil- 
ity the network exhibits |j,|^|. However, the Laplacian 
eigenvalues provide little insights into how the network 
topology affects the dynamics of synchronization. Thus, 
it is important and useful to investigate the connection 
between the synchronizability of networks and some in- 
tuitive graph invariants. Since the degree sequence plays 
a central role in revealing the scale-free property of net- 
works, its link with synchronizability has recently been 
studied |2]- It was found that the degree sequence itself 
is unable to determine the synchronizability of networks, 
although it does provide efficient information under cer- 
tain homogeneity constraints 8j. The diameter, which is 
defined to be the maximal length of the shortest paths be- 
tween any two vertices, has straightforward implication 
on the structure of a graph. Some previous experimental 
studies suggested that a small diameter might be advan- 
tageous to synchronization. This seems quite reasonable 



at the first sight, since a small diameter reduces the dis- 
tance of signal transmission on the network, but lacks 
theoretical verification. 

The purpose of this paper is to clarify the connection 
between the diameter and the synchronizability of net- 
works. Instead of using experimental methods directly, 
we employ both constructions and analytic estimations. 
We first explicitly construct a class of determinate net- 
works, named polarized networks^ which keeps the diam- 
eter invariant while the network size increases. Inequal- 
ities in spectral graph theory are used to show that A2 
and Aa/An approach zero as the network size increases to 
infinity, and thus these networks are poorly synchroniz- 
able. Then, we use Erdos-Renyi random graph model [9j 
to obtain a class of random networks with the same di- 
ameter, whose synchronizability proves to be quite good. 
These facts are somewhat surprising, since they point 
out that the diameter may not be an appropriate quan- 
tity to characterize the synchronizability of networks. To 
make this observation consistent with those experiments 
that suggested the advantage of a small diameter, we 
further explore the role that the diameter plays in the 
synchronization of networks. Analytic estimates on the 
synchronizability of networks in terms of the diameter are 
derived and show that the diameter does set lower and 
upper bounds for networks of a certain size. By these es- 
timates, we find that the synchronizability becomes more 
ficxible as the network size increases. We also present nu- 
merical results that support the theoretical analysis. 



II. NETWORK CONSTRUCTIONS 

In this section, we construct two classes of networks 
with prescribed diameter but very different synchroniz- 
ability. To do this, we first recall some concepts in graph 
theory [lj| . A complete graph is a graph with no loops in 
which each pair of vertices are joined by one edge. The 
complete graph of order n is denoted by Kn- Let A and 
B be two subsets of V{G), the vertex set of graph G, and 
define the distance between A and B to be the minimal 
length of all paths between a vertex in A and a vertex in 
B. 





FIG. 1: An example of polarized networks: P(26, 5). Two 
complete graphs K12 are joined by a path of length 3. 



A. Networks with Poor Synchronizability 

Suppose the network size n and the diameter D are 
given. We construct the polarized networks P{n, D) by 
joining two complete graphs K\ ^-n+a 1 and Kr n-p+s i by 

a path of length D — 2; an example is shown in Fig. 1. 
It is easy to see that P{n, D) has size n and diameter D. 
Note that the orders of the two complete graphs differ by 
at most one, and they are highly clustering at the ends 
of a single path, from which the word "polarized" comes. 
To derive an upper bound on the second smallest 
Laplacian eigenvalue A2 of P{n,D), we need a result in 
[ll|. If A and B are two subsets of V{G) with distance p, 
and F is the set of edges with at least one end vertex not 
in A or B, then we have the lower bound on the order of 
F: 



\F\ > p'\2 



\A\\B\ 
\A\ + \B\ 



(1) 



Now let A and B be the vertex sets of the two complete 
graphs in the construction of P{n,D). Then clearly, we 
have \A\ = [ii^fi^J , \B\ =jB^f+3] , |F| = Z? - 2, and 
p = D — 2. It follows from 



m 



2 

that 



Ao < 



\F\i\A\ + \B\ 
p^\A\\B\ 



D 



4(n 



-D 



2) [ii^f±2J p 



.-g+3" 
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{D-2){n-D + 2){n - D + A)' 



< 



{D - 2){n - D + iy 
4 
[D - 2){n - D + i)' 



D even 



D odd 



(2) 



On the other hand, a lower bound on the largest Lapla- 
cian eigenvalue A„ of a graph was given in [i3|. Let A 
be the maximal degree of a graph of order n. Then 



Now in P{n, D) 



A„, > 



A„> 



1 



L» + 3 



-A. 

Therefore, 
n{n — D 



> 



3) 



2{n-l) 



Combining Q and ||2J) gives 



< 



(71-1) 



(D - 2) n (n - Z? + 3)2 ■ 



(3) 



(4) 



Letting n ^ 00 in ((21 and © yields that both A2 and 
A2/A„ approach zero, which implies that the synchroniz- 
ability of polarized networks can be arbitrarily poor as 
the network size increases. What is more, we will show in 
the next section that the polarized networks are nearly 
the optimal class of networks with poorest synchroniz- 
ability among all the networks with the same network 
size and diameter. 



B. Netw^orks w^ith Good Synchronizability 

Instead of explicitly constructing a determinate net- 
work, we construct a class of random networks based on 
Erdos-Renyi random graph model. By this construction, 
we show that although there exist some extremal cases 
like the polarized networks, a general class of networks 
has fairly good synchronizability. 

A random graph G{n,p) consists of n vertices and the 
edges chosen independently with probability p. For the 
purpose of our construction, we need to fix the diameter 
of the graphs for all network size n. However, the diam- 
eter of G{n,p) may vary when n increases, depending on 
how we choose p as a function of n. T hus, we need a re- 
sult concerning the diameter of G{n,p) n^. Suppose the 
functions D = D{n) > 3 and Q < p — p{n) < 1 satisfy 
that logn/D — 3 log log n -^ 00, p^n^^^ — 21ogn -^ 00, 
and p^^^n^^^ — 21ogn -^ —00. Then G{n,p) almost 
surely, i.e. with probability tending to one as n tends to 
infinity, has diameter D. Here, let I? > 3 be a constant 
and 



P = 



(anlogn)i/(^^'^) 



(5) 



where < £ < 1 and a can be any positive number. It 
is easy to verify that all the conditions mentioned above 
are satisfied, which yields that a random graph G{n,p) 
with p defined in ((SJ almost surely has diameter D. 

The random graphs have a very narrow degree dis- 
tribution. In fact, we have the proposition |l4| . say- 
ing that ii p = p{n) satisfies that pn/ log n ^ 00 and 
(1 —p)n\ogn -^ 00, then almost surely all the degrees 
of G{n,p) are equal to (1 + o{l))np. It is clear that the 
assumptions are satisfied under (|SJ). Hence, almost surely 
all the degrees of G{n,p) with (jsj holding are near np. 

Now we are concerned with the Laplacian eigenvalues 
associated with G{n,p). The Laplacian L and the ad- 
jacency matrix ^ of a graph are linked by the relation 
L = K — A, where K is the diagonal matrix with diagonal 
entries being the vertex degrees. Due to the narrow dis- 
tribution of degrees in G{n,p), we have that L ~ npI—A, 
and so 



A,, 



np - 



A,- 



l,n. 



(6) 



where A^ denote the eigenvalues of A. Therefore, what 
we need to know is the distribution of eigenvalues of the 



adjacency matrix A. We know that this distribution fol- 
lows a semicircle law |l5l | , and a further result concerning 
all the eigenvalues Xi except Ai [la| says that 



max Ai(G'(n,p)) — 0{^/rvp). 

2<i<n 



(7) 



And also notice that Ai stands far away from the bulk of 
the spectrum, at np. Combining © and yields 

Xi ~ np + 0{y^), i = 2,n. 

With (O , we obtain the approximations for the two quan- 
tities characterizing the synchronizability 



A2~np= {anlogny^'^°-'\ 



(8) 



We can solve this inequality and obtain an upper bound 
on A2/A„ 



cosh 



< 



f coslr\n--l)_\ 
\ D-1 ) 



1 



Xo 



^n cosh( ^°'^^:'^r'M + i 



/ cosh-i(ra-l) \ 



(11) 



Now note that a simple upper bound on A„ was given 
El by 



Xn <n 



(12) 



with equality holding if and only if the complement of 
the graph is disconnected. Combining (|10|l . (|ll|l . and 
(|12|l gives the estimates as follows. 

Proposition: Given the network size n and the diam- 
eter D, the following estimates on X2 and A2/A„ hold: 



A2_ ^ np - 0{y/rvp) 
Xn ~ np + 0{y/np) ' 



(9) 



which implies that A2 ^ 00 and A2/A„ —i- 1 as n —>■ 00. 
By this we have shown that the class of random networks 
we constructed keeps the diameter invariant but exhibits 
rather good synchronizability. 



III. ANALYTIC ESTIMATIONS 

By the constructions given in the previous section, we 
see that the synchronizability of networks with the same 
diameter may differ significantly. They also suggest that 
the diameter should not be an appropriate quantity to 
characterize the synchronizability of networks, although 
this is somewhat conflictive with intuition. Naturally, 
one could ask whether or not there exists some relation 
between the diameter and the synchronizability. The an- 
swer is yes; we will now confirm this relation by deriving 
some estimates on the quantities A2 and A2/A„. These 
estimates are based on a series of inequalities in spectral 
graph theory. 

A lower bound on A2 in terms of the network size n 
and the diameter D \V^ pointed out that 



A2 > 



uD' 



(10) 



Chung [li 
bound 



estimated the diameter D by giving the upper 



D < 



cosh ^{n— 1) 



cosh 



-1 t \„ + \2 



A„-As 



which implies that 



D-K 



cosh (n — 1) 

u-l /l + A,/A„ 



cosh I 



< A2 < n ^ ^ 

'"'^ cosh 



1 



4 A2 
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n'^D - Xr, ~ 
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/cosh ^{n— 1) 
/ cosh" (n — 1) 



cosh 



{ D-1 ; + 



1 



(13) 



(14) 



Proof: The upper bound on A2 follows from (|ll(l and 
(|12|l . and the lower bound on A2/A„ follows from H1U|) 
and |(T3l. □ 

To gain more insights into the connection between the 
diameter and the synchronizability, we find an asymp- 
totic order of the upper bound on A2/A„. It is easy to 
calculate that 



cosh 



/cosh"i(n-l)\ -, 

( sV^j-l (2n)i/(c-i)-2 



cosh 



/cosh -^(n— 1) 

V S3I 



+ 1 



(2n)i/(c-i)+2' 



With the estimates and asymptotic order, it is easily 
seen that a larger diameter D implies both smaller lower 
bounds and smaller upper bounds on A2 and A2/A„. This 
observation coincides with the intuition that a larger di- 
ameter is likely to be harmful to synchronization. From 
this point of view, the assertion that a small diameter 
benefits synchronization could be true, although it is not 
true for individual networks. 

Another observation worthwhile to point out is that a 
larger network size n implies smaller lower bounds and 
larger upper bounds on both A2 and X2/Xn- In other 
words, the synchronizability of networks becomes more 
flexible as the network size increases. This fact suggests 
that the indicator role of the diameter in synchronization 
will substantially weaken in large-size networks. 

Let us return to the polarized networks for which the 
estimates on A2 and A2/A„ were given in Q and Q), 
respectively. The asymptotic orders of the upper bounds 
are A/{nD) and 8/{n'^D), respectively, for large n and 
D, and n ^ D. Thus, by the Proposition, we see that 
the polarized networks are almost the optimal class of 
networks with poorest synchronizability. 
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FIG. 2: The relationships between the network size n and 
the two quantities that characterize the synchronizabihty of 
polarized networks P{n,D). The diameter D is set to 5 all 
the time. 



FIG. 3: The relationships between the diameter D and the 
two quantities that characterize the synchronizabihty of po- 
larized networks P{n, D). The network size n is set to 200 all 
the time. 



IV. NUMERICAL RESULTS 

In this section, we carry out numerical experiments 
to explore the connection between the diameter and the 
synchronizabihty for the two classes of networks we con- 
structed. Note that the polarized networks are deter- 
minate while the random networks are, of course, not. 
So the former case is much simpler, while some average 
methods should be used to eliminate random perturba- 
tions in the latter case. 

In the polarized networks, all we need to determine 
before generating the network topology are the network 
size n and the diameter D. To investigate the influence of 
network size n and diameter D on the synchronizabihty of 
polarized networks, we assigned a set of values to the two 
parameters, and plotted the trends of the two quantities 
that characterize the synchronizabihty in Figs. 2 and 3. 
According to the lower bound estimates in (|13|l and (|14|l , 
combined with the upper bound estimates in |(5J and l@J, 
the relationships between the pairs of parameters we have 
shown in Figs. 2 and 3 are expected to be all linear. The 
numerical results are quite consistent with and support 
this observation. They confirm that the synchronizabihty 
of polarized networks is very poor, and decreases toward 
zero at least linearly as the networks size and diameter 
increase. 

Now we turn to the random networks G{n,p) with p 
defined in Q. It is worth noticing that although any 
< £ < 1 and a > ensure the convergence of diameter 
to the specified value, a careful selection of their values 
will accelerate the convergence so that networks with a 
particular diameter may be obtained for relatively small 
network sizes. In our experiments, we always set e = 0.5, 
but choose different a values for different diameters D. 
The results of experiments for D — A and a — ^ are il- 
lustrated in Fig. 4. In the experiments, ten random net- 
works were generated for each network size n, and the 
averages of A2 and \2l\n were computed respectively. 
The diameter converged fast to and stayed stably at the 
specified value. This is shown in the bottom of Fig. 4, 
where N denotes the number of networks achieving diam- 
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FIG. 4: The performance of random networks with p defined 
in 0. For the set of fixed parameters D = 4, e = 0.5, a = 5, 
ten random networks were generated for each network size n. 
The averages of the two quantities related to the synchroniz- 
abihty are shown in the top two graphs, while A*', the number 
of networks achieving diameter D out of the ten for each net- 
work size n is shown in the bottom. 



eter D out of the ten for each network size n. Meanwhile, 
the two quantities A2 and A2/A„ related to the synchro- 
nizabihty went up slowly. These trends are in agreement 
with the analytic estimations. However, it is difficult to 
obtain the accurate growing speeds of the two quantities 
by numerical experiments, since the required computa- 
tions increase very rapidly along with the network size. 
More accurate numerical results are shown in Fig. 5, 
where 50 random networks were generated for each set 
of parameters and averages of A2 and A2/A„ were com- 
puted respectively, so that the effects of random pertur- 
bations could be reduced remarkably. First, it is clear 
that the increasing diameter has pulled down the syn- 
chronizabihty of networks. This is consistent with the 
intuitive observation that a large diameter would impede 
communications between nodes and synchronization in a 
certain class of random networks, although the observa- 
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FIG. 5: The average performance of the class of random net- 
works we constructed, under different diameters and network 
sizes. Here e — 0.5, and a is set appropriately to ensure fast 
convergence to the specified diameter. For each set of param- 
eters, 50 random networks were generated and averages of the 
two quantities related to the synchronizability were computed 
respectively. 



tion is absolutely not true for determinate networks, for 
instance, a polarized network and a typical sample in the 
random networks with the same diameter, which have 
been shown to perform distinctly in synchronization. 

Another observation which can be seen from Fig. 5 is 
that the class of random networks has an obvious trend 
in the improvement of synchronizability as the network 
size increases, which turns out to be the same as what has 
been confirmed by Fig. 4. However, it is worth noticing 
that this trend becomes quite faint for a large diameter. 
This can be seen from the asymptotic orders in (jSJ and 
0. Restricted by the complexity of computation, we 
were only able to conduct the experiments in a relatively 
small range of network size n, but the increase of network 
size was pulled down considerably by a large diameter D, 
according to ||SJ) and Q. 



networks, we found that the networks with the same di- 
ameter may have very different synchronizability. Thus, 
the diameter itself is inappropriate to be a quantity that 
characterizes the synchronizability of networks. We also 
derived analytic estimates on the range of synchroniz- 
ability of networks with specified diameter and network 
size, which revealed that a larger size of network admits 
of a more flexible synchronizability. 

An application of the network constructions is the de- 
sign of network topology in communications, multi-agent 
systems, and other engineering problems. Since the po- 
larized networks are almost the optimal class of networks 
with poorest synchronizability, they are especially suit- 
able to be used in the design of easily desynchronizable 
networks. Moreover, the nature of explicit constructions 
makes it quite easy to apply. The random networks with 
fixed diameter are also useful in synchronization-related 
applications that have special requirements on the diam- 
eter of networks. The fast convergence to the specified 
diameter enables us to generate desired networks of rel- 
atively small sizes. 

A noteworthy endeavor in the literature is to explicitly 
construct networks with best synchronizability. Recently, 
methods in numerical optimization such as simulated an- 
nealing have been utilized to obtain nearly optimal net- 
work topology [23. It was suggested that the optimal 
network topology should possess several important prop- 
erties, including homogeneous degree, betweenness, and 
distance distributions, large girths, small diameters, no 
community structure, etc. However, the procedure to 
generate such a network topology remains an open ques- 
tion. As pointed out in this paper and previous research, 
individual properties such as small diameters, homoge- 
neous degree distributions, etc., are insufficient to op- 
timize the synchronizability of networks. A successful 
constructing algorithm should instead take into account 
a full combination of these favorable properties. 



V. CONCLUSIONS 

We investigated the connection between the diameter 
and the synchronizability of networks. By constructing 
determinate polarized networks and a class of random 
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